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Abstract 

We investigate strong pairing fluctuations and effects of a harmonic trap in the superfluid phase 
of an ultracold Fermi gas. Including amplitude and phase fluctuations of the inhomogeneous 
superfluid order parameter A(r) in a trap within a combined T-matrix theory with the local density 
approximation, we examine local properties of single-particle excitations and a thermodynamic 
quantity in the BCS-BEC crossover region. Below the superfluid phase transition temperature T c , 
we show that inhomogeneous pairing fluctuations lead to a shell structure of the gas cloud in which 
the spatial region where the ordinary BCS-type superfluid density of states appears is surrounded 
by the region where the pseudogap associated with strong pairing fluctuations dominates single- 
particle excitations. The former spatial region enlarges to eventually cover the whole gas cloud 
far below T c . We also examine how this shell structure affects the photoemission spectrum, as 
well as the local pressure. Since a cold Fermi gas is always trapped in a harmonic potential, our 
results would be useful for the study of strong-coupling superfluid physics, including this realistic 
situation. 
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I. INTRODUCTION 



Since the background physics of superfluid Fermi gases IrM] is similar to that of metallic 
superconductivity, the former system is now expected as a useful quantum simulator for 
the latter. A tunable pairing interaction associated with a Feshbach resonance {5 -10] in a 
cold Fermi gas enables us to study Fermi superffuids from the weak-coupling BCS (Bardeen- 
Cooper-Schrieffer) regime to the strong-coupling BEC (Bose-Einstein condensation) regime 
in a unified manner [5Nl5|. In this BCS-BEC crossover, one can systematically examine 
strong-coupling effects by adjusting the interaction strength. 

While there exist various similarities between a superfluid Fermi gas and metallic su- 
perconductivity, the presence of a trap potential is peculiar to the former. Because of this 
confined geometry, physical quantities naturally become inhomogeneous. Thus, when an 
experiment has no spatial resolution, it gives spatially averaged data. For example, the 



photoemission-type experiment developed by JILA group [16l4l9j] so far has no spatial reso- 
lution, so that the observed data involve single-particle excitation spectra at various spatial 
positions. Thus, the spatial inhomogeneity is a crucial key in considering single-particle 
properties of a superfluid Fermi gas by using this experiment. Since a bulk superconductor 
is usually a uniform system, this problem is also important in using a superfluid Fermi gas 
clS db quantum simulator for superconductivity. 

In this paper, we investigate effects of a harmonic trap on strong-coupling properties of 
a superfluid Fermi gas in the BCS-BEC crossover region. In particular, as a typical strong- 
coupling phenomenon, we deal with the pseudogap problem 20M27I] . In this phenomenon, 
strong pairing fluctuations induce a dip structure in the density of states above the superfluid 
phase transition temperature T c . In our previous paper for a uniform Fermi gas 28], we 
showed that the pseudogap still remains just below T c . Then, in a trapped superfluid Fermi 
gas, we can expect the inhomogeneous situation that while the BCS-type superfluid density 
of states appears in the trap center, the pseudogap is still dominant in the outer region of 
the gas cloud (where the superfluid order parameter is small and pairing fluctuations are 
strong) . 

To confirm such a shell structure, it is convenient to examine the superfluid local density of 
states (LDOS), as well as the local spectral weight (LSW). In this paper, using a combined 
T- matrix theory with the local density approximation (LDA), we identify the pseudogap 
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region in the phase diagram with respect to the interaction strength, temperature, and 
spatial position. We also discuss how the shell structure affects the photoemission spectrum. 
Besides single-particle excitations, thermodynamic properties has also attracted much 



attention in the BCS-BEC crossover region [8, 



29N36j| . In this paper, we consider the local 



pressure P{r) as a typical thermodynamic quantity. Recently, Ho and Zhou[34| proposed 



a useful idea to determine P(r) from the density profile. Using this, Nascimbene and co- 
workers}^ measured P{r) in the unitarity limit of a 6 Li Fermi g function of the 
LDA fugacity ((r) = e^' kl}T (where /i(r) is the LDA chemical potential). They reported 
that the observed pressure is well described by the Fermi liquid theory. In this paper, we 
clarify to what extent strong pairing fluctuations affect P(r), comparing our results with 



the experimental data 



36]. 



This paper is organized as follows. In Sec. II, we explain our formulation. In Sec. Ill, we 
present our numerical results on LDOS and LSW in the BCS-BEC crossover below T c . We 
also identify the region where the pseudogap dominates single-particle excitations in the 
phase diagram of a trapped Fermi gas. In Sec. IV, we consider the photoemission spectrum. 
We clarify how the inhomogeneous pseudogap affects this quantity below T c . In Sec.V, we 
treat the local pressure P(r). We show that the calculated P(r) agrees well with the recent 
experiment on a 6 Li Fermi gas 36|. Throughout this paper, we set h = &b = 1. 



II. FORMULATION 

We consider a two-component superfluid Fermi gas in a harmonic pot ential, described by 
the BCS Hamiltonian. In the Nambu representation, it has the form 28l. 1371]. 



H 



En[«, 



t 3 - An 



%-U^2p+(q)p-(-q). 



(Effects of a harmonic trap is later included within LDA.) Here, 



(1) 



c 



Cpt 

-pi 



(2) 



is the two-component Nambu field, where <3 p a is the creation operator of a Fermi atom with 
pseudospin a (=t, \), describing two atomic hyperfine states. £ p = e p — p = p 2 /(2m) — p 
is the kinetic energy, measured from the Fermi chemical potential \x (where m is an atomic 
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FIG. 1: (a) Self-energy correction Y, p (iu) n , r) in the combined T-matrix theory with LDA. The wavy 
line is the particle-particle scattering matrix T q ' s (iv n ,r) (s,s' = ±), describing amplitude and 
phase fluctuations of the superfluid order parameter. The solid and dashed lines describe the LDA 
mean-field Green's function Gp(iu n ,r) = [ico n — £ p (r)r3 + A(r)Ti] _1 , and the pairing interaction 
—U, respectively. The bubble diagram describes the pair-correlation function Hq S (iv n ,r). 

mass). —U (< 0) is a pairing interaction, which is assumed to be tunable by adjusting the 
threshold energy of a Feshbach resonance. Tj (j = 1,2,3) are the Pauli matrices, acting on 
the particle-hole space. In this paper, we take the superfluid order parameter, 

A = uJ2(c- P i^), (3) 
p 

to be real and be parallel to the Ti-component. In this case, the generalized density operators 
PM) = E P ^p+q/2 r i^p-'?/2 (J = 1,2) in p±(q) = [p^q) ± ip 2 (g)]/2_describe amplitude 
and phase fluctuations of the superfluid order parameter, respectively ^ 



As usual, we measure the interaction strength in terms of the s-wave scattering length 



a s , which is related to the pairing interaction —U as 15] 



where u c is a high-energy cutoff. 

We now include effects of a harmonic trap potential V(r) = mu^ r r 2 /2 within LDA. 
This extension is achieved by simply replacing the Fermi chemical potential p by the LDA 



expression p{r) — fi — V(r)[2Q |39(. Various quantities then depend on the spatial position r 



through fi(r). For example, the 2 x 2-matrix LDA single-particle thermal Green's function 



4 



is given by 



G p (iu n ,r) = - — — — — — - 

iu n - £p{r)T 3 + A{r)Ti - 2J p (iu n ,r) 



(5) 



where £ p (r) = e p — fi(r), and u n is the fermion Matsubara frequency. A(r) is the LDA super- 
fluid order parameter. The LDA self-energy T, p (iu n , r) involves effects of pairing fluctuations 
within the T- matrix approximation, which is diagrammatically given by Fig. H|22L |28 | . Sum- 
ming up these diagrams, we obtain 



S p (zo; n ,r) = ~T^2 ^2 T q S '( iu n,r)T_ s G p+q (iui n + iu n ,r)T_ s/ . 

q,v n s,s'=± 



(6) 



zu, r 



Here, r± = [ti ± ir 2 ]/2, and v n is the boson Matsubara frequency. G p { 
£ P ( r ) r 3 + A(r)ri] _1 is the LDA mean-field Green's function. The particle-particle scattering 
matrix r q ' s '(iu n , r) has the form 



F+ (zz/ n ,r) T+ + (iu n ,r) 
r„~On,r) T- + (iu n ,r) 



-U 



i -i 



, , ( n q + (zz/ n ,r) Il++(iu n ,r) 
II~(ii/ n ,r) n+~(zz/ n ,r) 



(7) 



where 



(8) 



is the lowest-order pair correlation function with respect to the pairing interaction. Carrying 
out the c<; n -summation in Eq. (JSJ), we obtain 

sA(r) 2 E p+q/2 (r) + sE p _ q/2 {r) 



4 s t^i ^ E p+q/2 (r)E p - q/2 (r) z/ 2 + (£ P+<7 / 2 (r) + sE p . q/2 (r)) 2 



x 



tanh 



p+q/2l 



2T 



+ s tanh 



2T 



(9) 



> = zEE 



S=±l p 



1 + S 



+ 



1 - 



X 



Zp+q/2(r)Zp- q /2{r) 

E p+q/2 (r)E p _ q/2 (r) 
E p - q / 2 {r) 



l-s 



E p + q / 2 (r) 

tanh^^ 2(r) U g tanh^^ 2(r) 



iu n - (E p+q/2 (r) + sE p _ q/2 (r)) 

iv. 

n 



vt + (E p+q/2 (r) + sE p _ q/2 (r)y 



(10) 



\ 2T j \ 2T 

The other components are given by U q {iu n ,r) = n+ + (iz/ n ,r), and IT~ + (2z/ n ,r) = 
. In Eqs. © and (JTOP, E p (r) = ^(r) 2 + A(r) 2 is the LDA Bogoliubov 



n+- 



-iv n ,r 



single-particle excitation spectrum. 
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FIG. 2: (Color online) (a) Calculated LDA superfluid order parameter A(r = 0). The inset shows 
T c , normalized by the Fermi temperature Tp. (b) Spatial variation of A(r). Rp = y / 2ep / (mco^) 
is the Thomas Fermi radius, where £f is the Fermi energy, (c) Chemical potential \x. The dashed 
line shows T c . The inset shows fi(T = T c ). 

We self-consistently determine the local superfluid order parameter A(r) and the chemical 
potential /i, by solving the LDA gap equation, 

1 = UY — L-tanh^P, (11) 
^ 2E p (r) 2T ' v ' 

together with the equation for the total number of Fermi atoms, 

/■OO 

N = 4irr 2 drn(r), (12) 
Jo 

where n{r) = 2T J2 P ui n ^ p (iuj n , r) |ne* 5tJ " is the particle density. Within the framework of 
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(a)(W ! = 



-1 




FIG. 3: (Color online) Atomic density profile n(r) in the BCS-BEC crossover below T c . (a) 
(kpas)' 1 = -1 (BCS side), (b) (/c F as) _1 = (unitarity limit), (c) (/cpOs)" 1 = 0.8 (BEC side). 



LDA, T c is determined from the BCS-type T c -equation in the trap center (r = 0) 25 



26 



l = f/V — tanh^. (13) 
p y 

Figure [2] shows the self-consistent solutions of A(r) and //. We will use these numerical 
results in calculating various quantities in later sections. 

We note that the LDA gap equation (jlip gives the vanishing superfluid order parameter 
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FIG. 4: (Color online) Mean-filed contribution riMF( r ) and fluctuation contribution n^(r) to the 
density profile n(r) at 0.5T C in the unitarity limit. The LDA superfluid order parameter A(ro) 
vanishes at r = ro, at which a cusp appears in ns(?"), . 

A(r) = for r > r , where r is determined from the equation, 

1 = C/V — L_tanh^M. (14) 
^2£ p (r ) 2T 1 } 

However, this vanishing A(r > ro) is an artifact of LDA, because the superfluid order 
parameter should be finite everywhere in the gas cloud below T c . Since the correct A(r) 
would be large around the trap center, r obtained in LDA should be interpreted as a 
characteristic radius, inside of which the magnitude of A(r) is large. 

As pointed out in Ref. 2l| . the LDA atomic density profile n(r) has a cusp at r = ro, which 
becomes more remarkable as one approaches the strong-coupling regime, as shown in Fig. 
[3j Dividing n(r) into the sum of the mean-field part «mf('") = 2TJ2 puJn Gp(w n , r)\ue l5uJn 
and the fluctuation contribution, 

ra B (r) = [G P (co n ,r) - G° p (u) n ,r)] \ n e i5 ^. (15) 

we find in Fig. H]that the cusp only appears in riR fr) . In a Bose gas BEC, the LDA density 
profile at T c is known to exhibit a cusp at r = 0[40], which is also seen in Figs. EJ^b) and 
(c) at T c . Since ne(r) reduces to twice the molecular density profile in the BEC limit, the 
cusp singularity seen in Fig. [3] is found to be the same artifact of LDA as in the Bose gas 
case j4o| . 

Once A(r) and fi are determined, the LDA single-particle spectral weight (LSW) 
A(p,u,r) is conveniently calculated from the analytic continued Green's function as 

A(p, cu, r) = lmG p (iu n —> u + 18, r)|n. (16) 

71 



The local density of states (LDOS) p(u, r) is related to LSW as 

p(u,r) = ^A(p,uj,r). (17) 
p 

For the infinitesimally small positive number S appearing in Eq. [161 to avoid unphysical 
fine structures in LDOS and LSW, we take 5 = O-Olew in numerical calculations. 

The .oca. spectra! weight (LSW) i, a ls0 r e.ated to the photoe misS io„ spectrunfi Q 
In this experiment, atoms in one of the two atomic hyperfine states (= | f)) are transferred 
to another hyperfine state |3) (^ | f), | D) by rf-pulse, and one measures the rf-tunneling 
current from | t) to 13). In 40 K Fermi gases, because the so-called final state interaction 
can be safely ignored jig. \v\, |3) may be treated as a non-interacting state. Using the linear 
response theory, we obtain the rf-tunneling current as, 

I(p, Q, r) = 2irt 2 F A p (£ p (r) - Q, r)/(£ p (r) - Q). (18) 



For the derivation of Eq. f[T8l . we refer to Ref. 26]. tp represents a transfer matrix element 
between | f) and |3). f(e) is the Fermi distribution function. Noting that the current 
photoemission-type experiment does not have spatial resolution, we take the spatial average 
of Eq. (TUD as 

WP,0) = J drA p (Ur) ~ n,r)/&(r) - «)• (19) 

Here, Rp = \/2ef j (mu^ v ) is the Thomas- Fermi radius (where Ef is the Fermi energy), and 
Q = u^—us is the energy difference between the incident photon energy cul and the en ergy cu.s 



of the final state |3). Equation (1T9"|) is related to the observed photoemission spectrum [lq.ll 7] 
p 2 A p {oj) f (u) , as well as the occupied density of states p(o;)/(o;), as 



p 2 A p (uo)f(u) = p 2 h vc (p, n^t P -ou), (20) 
1 

2~' 2 



= — / dpp 2 A p {u)f{u). (21) 



To calculate the local pressure in a gas cloud, we employ the idea proposed by Ho 
and Zhou34(. That is, using the Gibbs-Duham equation for the local pressure, dP(r) = 



n(r)/j,(r) + s(r)dT (where s(r) is the entropy density), we obtain dP(r) = n(r)dfi(r) for a 



fixed value of T. Then, using the relation dfi(r) = mu 2 rdr and P(r — > oo) = 0, one finds 34j 

P(r) = mu 2 [ r'dr'n{r'). (22) 



Substituting the calculated particle density n{r) shown in Fig. [3] into Eq. (|22|) . we numeri- 
cally carry out the integration in Eq. (T22]) to obtain the local pressure P{r). 
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FIG. 5: (Color online) Calculated local density of states p(oj, r) in the unitarity limit ((k F a s ) 1 = 0). 
(al)-(a4) r = 0. (bl)-(b4) r = 0.5R F . (cl)-(c4) r = R F . 

III. PSEUDOGAP AND SUPERFLUID GAP IN A TRAPPED FERMI SUPER- 
FLUID 

Figure [5] shows the local density of states (LDOS) p(u, r) in the unitarity limit of a 
superfluid Fermi gas. In the trap center (panels (al)-(a4)), a large dip structure associated 
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with strong-pairing fluctuations is seen around uj = at T c . This pseudogap has already 



appeared above T c [26[ (although we do not show the result here). At 0.9T C (panel (a2)), 
although the LDA superfluid order parameter A(r = 0) is finite, LDOS still has a finite value 
at uj — 0, because of residual pairing fluctuations below T c . The BCS-type full gap structure 
can be only seen below T ~ 0.8T C . Far below T c (panel (a4)), the ordinary BCS-type density 
of states is realized, being accompanied by sharp coherence peaks at the excitation edges. 
Panels (al)-(a4) indicate that the pseudogapped LDOS at r = smoothly changes into the 
superfluid density of states with decreasing the temperature. 

At r = 0.5-Rp; since pairing fluctuations are weaker than those in the trap center, the 
pseudogap (dip) structure does not appear at T c , as shown in Fig. EJ^bl). As mentioned 
previously, the LDA superfluid order parameter A(r) vanishes at r > r (T), so that pairing 
fluctuations at r > r continues to develop even below T c . Because of this, the gradual 
development of the pseudogap around uj = is seen in Fig. Mjo2) and (b3). (We note that 
A(r = 0.51-Rf) = at T = 0.7T C .) Far below T c , since the LDA superfluid order parameter 
at r = 0.5i?F becomes finite, the ordinary BCS superfluid density of states is obtained, as 
shown in Fig. |5]^b4). 

The particle density is very low around the edge of the gas cloud (r ~ R F ), so that the 
superfluid order parameter is small and pairing fluctuations are weak there. As a result, 
LDOS at r = R F shown in panels (cl)-(c4) is close to the density of states of a free Fermi 

gas, 



P {uj, r = R F )~ v / u + n( r )G(u + /x(r)), (23) 

where Q(x) is the step function. 

Since pairing fluctuations are weak in the BCS regime, the ordinary BCS-type superfluid 
density of states is soon realized below T c . (See Figs. |6]^al)-(a4).) On the other hand, LDOS 
in the BEC regime (Figs. [6]^bl)-(b4)) is characterized by a fully gapped structure, reflecting 
the large binding energy of a tightly bound molecule, which has been already formed above 
T c . While the gap size is almost unchanged below T c , the growth of the coherence peaks can 
be seen with decreasing the temperature. 

for a uniform Fermi gas, we have introduced two characteristic 



In our previous paper 



temperatures T* and T* to identify the pseudogap regime. T* corresponds to the ordinary 
pseudogap temperature, namely, it is defined as the temperature at which a dip (pseudogap) 
structure appears in the density of states above T c . At T*, the superfluid density of states 
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FIG. 6: (Color online) Calculated local density of state at r = 0. (al)-(a4) (kpa s ) 1 = — 1 (BCS 
regime). (bl)-(b4) (kpas)' 1 = 0.8 (BEC regime). 

p(oj = 0) is suppressed by 50% compared to the value of p(ui = 0) at T c . Since the mean- 
field BCS state always has the vanishing density of states at u> = below T c , the fact that 
p(u = 0) still has a large value at T* < T < T c means the importance of pairing fluctuations 
there. We have regarded the region T* < T < T* as the pseudogap (PG) regime [281 . where 
pairing fluctuations dominate single-particle excitations. 
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FIG. 7: (Color online) (a) Three-dimensional phase diagram of a trapped superfluid Fermi gas. 
The solid line at r = is T c . For the definitions of T*(r) and T*(r), see the text. ro(T) is the 
spatial position where the LDA superfluid order parameter vanishes at T. 'SF' is the superfluid 
region where single-particle excitations are close to the ordinary BCS-type. 'PG' is the pseudogap 
regime characterized by the pseudogapped local density of states. 'NF' is the normal Fermi gas 
region where neither the superfluid gap nor the pseudogap appears in LDOS. In the BEC regime 
where [i < 0, we also plot 2\/j,(T c )\, which gives the characteristic temperature where two-body 
bound states appear. The right side of this line may be viewed as a molecular Bose gas, rather 
than a Fermi gas. In this regime, 'MBEC' is the region which is well described by the BEC of 
tightly bound molecules. 'NB' is the region of a non-condensed molecular Bose gas. We also 
show the tomographic views of this phase diagram in panels (b) (/cfOs) -1 = —0.5 (BCS side), (c) 
(&FQs) 1 = (unitarity limit), and (d) (kpa^^ 1 = 0.5 (BEC side). 

Extending the above discussion to the present trapped case, we introduce two charac- 
teristic temperatures T*(r) and T*(r). T*(r) is denned as the temperature at which a dip 
(pseudogap) structure appears in p(u,r) around u — 0. T*(r) is determined from the con- 
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FIG. 8: (Color online) Spatial variation of LDOS at T = 0.5T C in the unitarity limit. 

dition that p(u> = 0,r) is suppressed by 50% compared with the value at the temperature 
where the LDA superfluid order parameter at r becomes finite. 

Using T*(r) and T*(r), we obtain the phase diagram of a trapped Fermi gas in Fig. 



re overall structure is essentially the same as the phase diagram for a 
. As in the uniform case, we call the region between T*(r) and T*(r) 



[3(a). At r = 0, t 
uniform Fermi gas 

the pseudogap (PG) region. The BCS-type superfluid density of states only appears below 
T*(r) (SF region). Above T*(r) (NF region), the pseudogap is absent in LDOS, where 
excitation properties are close to those of a normal Fermi gas. 

In Fig. Ufa), we also plot 2|/i(T c )| in the BEC regime where the Fermi chemical potential 
n(T = T c ) is negative. Since 2|//| reduces to the binding energy E\> = l/{ma 2 s ) of a two-body 
bound state in the BEC limit, this line physically describes the characteristic temperature 
where two-body bound molecules appear, overwhelming thermal dissociation. Thus, below 
T ~ 2\fi(T c )\, the system may be viewed as a molecular Bose gas, rather than a Fermi atom 
gas. In this strong-coupling regime, T c is well described by the BEC phase transition from 
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a normal-state Bose gas (NB) to the molecular BEC (MBEC). 

We briefly note that, in Fig. [7J T c is only the phase transition temperature. T*(r), T*(r), 
and 2\n(T c ) \ are crossover temperatures without being accompanied by any phase transition. 

Figures [7(b)- (d) show that the PG region (T*(r) < T < T*(r)) always exists along the 
r (T)-line determined by Eq. (|14p . Since the LDA order parameter A(r) becomes finite 
below this line, in a sense, the ro(T)-line may be interpreted as the "local superfluid phase 



transition temperature (= T c (r))" within the LDA picture 42|. Thus, these panels indicate 
that pairing fluctuations at r become strong near T c (r), leading to the pseudogap around 
ro(T). 

At a fixed temperature below T c , Figs. [3(b)- (d) indicate that a trapped Fermi gas exhibits 
a shell structure. For example, in the unitarity limit at T = 0.5T C , one finds from panel (c) 
that the SF, PG, and NF regions occupy the spatial regions, < r < 0.55-Rf, 0.55-Rf < r < 
0.67-Rp, and 0.67-Rp < r , respectively. In this shown in Fig. [SJ while the BCS-type 

LDOS is obtained when r < 0.55-Rf, the pseudogap is seen at r = r = 0.64-Rp. However, 
the pseudogap does not appear when r > 0.67-Rf- 

As mentioned previously, the vanishing A(r) for r > ro is an artifact of LDA. In this 
regard, we note that the PG region also exists below the ro(T)-lines in Figs. E^b)-(d), where 
A(r) is small but finite. Thus, the PG region is expected to exist, even when one includes the 
finite value of A(r > r ) by a more sophisticated treatment. In such an improved theory, 
since A(r > r ) suppresses pairing fluctuations to some extent, the PG region would be 
narrower than the LDA result. 

The "SF-PG-NF shell structure" can be also seen in the spatial variation of local spectral 
weight (LSW) A p (co,r) in the unitarity limit. In the trap center (Fig. [9(a)), the gapped 
spectral structure is close to the ordinary BCS-type spectral weight, 



A* p cs (u, r = 0) = J I (l + 5(co - E p ) + J I (l - 6(u + E p ), (24) 

where E p = \fi P + A 2 (0) is the Bogoliubov single-particle excitation spectrum at r = 0. 
Since A(r) is smaller in the outer region of the gas cloud, the gap size seen in LSW also 
becomes small, as shown in panel (b). However, although A(r) vanishes at r = r = 0.64i?F, 
A p ((jj,r) in panel (c) still has a gap-like structure in the low momentum region, which is 



characteristic of the pseudogap phenomenon 24], l25| . This pseudogap becomes obscure in 
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FIG. 9: (Color online) Intensity of local spectral weight (LSW) A p (u,r) at T = 0.5T C in the 
unitarity limit, (a) r = 0. (b) r = 0.6R F . (c) r = 0.64R F (= r ). (d) r = 0.7i? F . (e) r = R F . The 
intensity is normalized by sZ . The same normalization is also used in Fig. [TUJ 

panel (d), to eventually disappear in panel (e). At r = Rp, the spectral peak line is close to 
the free-particle dispersion uj = p 2 /(2m) — /x(i2 F ). 

In the weak-coupling BCS regime, while the BCS-type gapped spectrum, as well as the 
free-particle-like peak line, can be seen at r = and r = i? F , respectively, a pseudogap 
structure does not appear at r = r , as shown in Fig. [TW al)-(a3). In the uniform case, it 
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FIG. 10: (Color online) Intensity of LSW at T = 0.5T C . (al)-(a3) {hpa^' 1 = -1 (BCS regime). 
(bl)-(b3) (fcpas)- 1 = 0.8 (BEC regime). 



has been shown [24] that the pseudogap in the BCS regime can be more clearly seen in the 
density of states than the spectral weight, leading to different pseudogap temperatures that 
are determined from these quantities. In particular, when (/c F a s ) _1 = —1, the pseudogap is 
almost invisible in the spectral weight even at T C [24J, although a dip structure appears in 
the density of states. In the present trapped case, a similar situation is considered to occur 
in panel (a2), although the pseudogap region exists along the ro-line in the phase diagram 
in Fig. [7] (which is obtained from the local density of states). 

In the BEC regime with /i < 0, Figs. [TUf bl)-(b3) show that, while the lower peak line 
gradually disappears with increasing r, the upper branch reduces to u = p 2 /(2m) — /x(r) 
around the edge of the gap cloud. The upper branch is related to the dissociation of tightly 
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FIG. 11: (Color online) Calculated intensity of photoemission spectrum p 2 A p (uj) f (to) . (al)-(a3) 
(^a s ) _1 = —1 (BCS regime). (bl)-(b3) (^a s ) _1 = 0.8 (BEC regime). The intensity is normalized 
by 27rf§/(2m). The same normalization is used in Fie;. [12] 



bound molecules that have been already formed above T c , so that it exists even for A(r) = 0. 
On other hand, since the lower branch is associated with a particle-hole coupling induced 
by the superfluid order parameter 2J, |25(, it is absent when r > r (T) (where A(r) = 0). 
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FIG. 12: (Color online) Calculated intensity of photoemission spectrum p 2 A p (lu) f (lu) in the uni- 
tarity limit (/cfOs) -1 = 0. The solid line is the free-particle dispersion uj + fj, = p 2 /(2m). The 
dashed line and dotted line show peak positions of the spectrum. 
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IV. PHOTOEMISSION SPECTRUM IN A TRAPPED SUPERFLUID FERMI 
GAS 



Figures [TT1 and [T21 show photoemission spectra p 2 A p {uj) f (co) . As discussed in Ref. [251]. 
the calculated spectra at T c (Figs. ITTT al'l. ll(bl) and IT27 a)) agree well with the recent 



experiment on a 4U K Fermi gas[16|, |17|. That is, starting from the weak-coupling BCS 
regime, a sharp peak line along the free particle dispersion u + \x = p 2 / (2m) in Fig. ITTT al) 
becomes broad in the unitarity limit (Fig. IT27 al)). which eventually splits into an upper 
sharp branch and a lower broad branch in the BEC regime (Fig. ITTT bl)). 

In the weak-coupling BCS regime, the overall spectral structure almost remains un- 
changed below T c , as shown in Fig. [TTT al)-(a3). However, when we carefully look at panel 
(a2), we find the splitting of the spectral peak around p/k F = 1. The upper peak is along 
the free particle dispersion to + \i — p 2 /(2m), so that it comes from the NF region around 
the edge of the gas cloud. On the other hand, the lower peak line does not exist in panel 
(al), and the momentum dependence is similar to the hole branch of the BCS Bogoliubov 
excitation spectrum (which is given by E p = - v /(p 2 /(2m) - e F ) 2 + A 2 in the BCS theory). 
Thus, the lower peak is considered to come from the SF region around the trap center. In 
this sense, the spectral structure seen in Fig. ITTT a2) is consistent with the shell structure 
discussed in Fig. [7J Since the SF region spreads out to the whole the gas cloud far below 
T c , the upper peak line disappears in Fig. [TTT a3). 

In the unitarity limit at T c , while the peak line along the free particle dispersion u> + fi = 
p 2 / (2m) in Fig. [T27 a) comes from the NF region, the broad spectral structure below this 



reflects the PG region 



25 



26|. Since the SF region only appears below T*(r = 0) = 0.81T C 



(See the phase diagram in Fig. [71), the photoemission spectrum almost remains unchanged 
at T = 0.9T C , as shown in Fig. [T2T b). The SF region starts to develop from the trap 
center below T*(r = 0), so that the spectral weight gradually move to the lower peak line, 
as shown in panels (c)-(e). In this low temperature region, in addition to the hole branch 
of Bogoliubov excitations, one also slightly sees the peak line corresponding to the particle 
branch of Bogoliubov excitations below the free particle dispersion (dotted line in panels (c) 
and (d).) Far below T c , because the SF region covers the whole gas cloud, and because the 
Fermi distribution function in Eq. (1191) suppresses the spectral intensity in the high energy 
region, the photoemission spectrum is dominated by the lower Bogoliubov branch, as shown 
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FIG. 13: (Color online) Calculated occupied density of states p(uS)f(u) in the superfluid state below 
T c . (a) (kFds)- 1 = -1 (BCS regime), (b) (k F a a )- x = (unitarity limit), (c) (&fci s ) _1 = -0.8 
(BEC regime). 

in Fig. ff2Te). 

In the strong-coupling BEC regime, Figs. [TTT bl)-(b3) show that the photoemission spec- 
trum already splits into an upper and lower branches at T c , reflecting a large molecular 
binding energy. Since the upper branch is suppressed by the Fermi distribution function in 
Eq. ( II 9p far below T c , the lower branch is only seen in Fig. H~lT b3). 

Figure [131 shows the occupied density of states p(cu)f(uj) in a trapped superfluid Fermi gas. 
In the BCS regime, as expected from the weak temperature dependence of the photoemission 
spectrum in Figs. ITl7 al)-(a3). the occupied density of states p{u)f{u) in Fig. ITW a) almost 
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FIG. 14: (Color online) Calculated local pressure P{r) in the BCS-BEC crossover regime of a 
trapped Fermi gas. (a) (/cf^s) -1 = —1 (BCS regime), (b) (kpa s ) _1 = (unitarity limit), (c) 
(/cf^s) 1 = 1 (BEC regime) . In panel (d), we plot P(r) as a function of n(r)/T. Experimental 
results on a 6 Li Fermi gas[36| are also shown in panel (d) (solid circles). Pq = Yg^-(2m) 3/ ' 2 e F //2 is 
local pressure of a free Fermi gas at r = and T = 0. Po(r) = rncu 2 f^r'dr' ^2 f(^ p (r')) is the 
local pressure of a free Fermi gas where the same values of T and \x as those used in calculating 
P(r) are taken. 

remains unchanged below T c . In contrast, in the unitarity limit (panel (b)), a peak starts 
to grow at uj/ef ~ —1 below T < 0.7T C to become a sharp peak far below T c . From 
the comparison with Fig. [T2T c). we find that this peak corresponds to the hole branch of 
Bogoliubov single-particle excitations. The growth of this sharp peak can be also seen in 
the BEC regime, as shown in Fig. [T3l c) . 
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V. LOCAL PRESSURE IN A TRAPPED FERMI GAS 



Figures IT4T a)-(c) show the local pressure P(r) in the BCS-BEC crossover regime of a 
trapped superfluid Fermi gas. In the BCS regime, panel (a) shows that P(r) is almost T- 
independent. Since P(r) is related to the particle density n{r) as Eq. (I2"2l [43] . this result 
reflects the weak temperature dependence of n(r) in this region. (See Fig. [3]^a)). When 
the interaction strength becomes strong, Figs. [3](b) and (c) indicate that particles tend to 
cluster around the trap center below T c . This enhances the local pressure P(r) around r = 0, 
as well as the decrease of P(r) around r = Rp, below T c , as shown in Fig. [T4T b) and (c). 

In the unitarity limit, since a s diverges, the system can be described by the single dimen- 
sionless parameter p(r)/T 44|. This university also holds in the present combined T-matrix 
theory with LDA, as discussed in the Appendix. Indeed, when we plot P(r) as a function 
of /i(r)/T, all the results at different temperatures are well fitted by a universal function, as 
shown in Fig. [T4T d). 

In panel (d), we also compare our result with the recent experiment on a 6 Li Fermi gas 
done by ENS group [36(. Without introducing any fitting parameter, our result is in good 
agreement with the experiment. 

We note that, although our theory correctly include the pseudogap effect associated with 
strong pairing fluctuations, one cannot see a clear signature of this strong-coupling effect in 



Fig. El Indeed, Ref. [36] reports that their experimental data can be well described by the 
Fermi liquid theory. From these, one finds that the local pressure P(r) is not sensitive to 
the pseudogap phenomenon. The LDA particle density n{r) in Eq. (|2"2"|) is given by 

/oo 
duf{u)p{u,r). (25) 
-oo 

Thus, even when the LDOS p(u, r) has a pseudogap structure around lo — 0, it would be 
smeared out to some extent by the ^-integration in Eq. ( 1251) . The pseudogap effect would be 
further smeared out by the spatial integration in Eq. (122]) . Because of these two integrations, 
the detailed pseudogap structure in the low energy density of states is considered to be not 
crucial for P{r). 
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VI. SUMMARY 



To summarize, we have discussed pseudogap phenomena and effects of a harmonic trap in 
the BCS-BEC crossover regime of a superfiuid Fermi gas. Extending our previous work for 
a uniform system to include effects of a harmonic trap within the local density approxima- 
tion (LDA), we clarified the coexistence of the superfiuid (SF) region where the BCS-type 
superfiuid density of states appears and the pseudogap (PG) region which is dominated by 
pseudogap phenomenon even below T c . From the spatial and temperature dependence of 
the local density of states, we identified the pseudogap region in the phase diagram with 
respect to the temperature, interaction strength, and spatial position. 

We have discussed the photoemission spectrum in a trapped superfiuid Fermi gas. In the 
BCS-BEC crossover region, the photoemission spectrum was shown to be strongly affected 
by the shell structure of a trapped superfiuid Fermi gas which consists of the SF region, PG 
region, and the normal Fermi gas (NF) region. Since the inhomogeneity and strong-pairing 
fluctuations are important key issues in considering a real trapped Fermi gas, our results 
would be useful for the study of strong-coupling effects of this system, including the realistic 
situation. 

We have also examined pseudogap effects on the local pressure P{r). While our strong- 
coupling result agrees well the recent experiment on a 6 Li Fermi gas done by ENS group 36], 
we showed that this thermodynamic quantity is not sensitive to the pseudogap appearing 



in the single-particle density of states. This is consistent with the statement [36J] that the 
observed pressure in the unitarity limit can be well described by the Fermi liquid theory. 
Since the pressure is not directly related to the detailed single-particle excitations compared 
with the photoemission spectrum, it is possible to occur that, while the pseudogap can be 
observed in the latter, such a strong-coupling phenomenon does not clearly appear in the 
former thermodynamic quantity. 

In this paper, we have treated the inhomogeneity of the system within LDA. While LDA 
has succeeded in explaining various properties of trapped Fermi gases, it cannot correctly 
describe the feature that the superfiuid order parameter is finite everywhere below T c . To 
overcome this, one needs a more sophisticated inhomogeneous strong-coupling theory than 
LDA. Since the presence of a trap potential is unique to the cold Fermi gas system, this 
problem would be an important challenge for the further development of the BCS-BEC 



24 



crossover theory in cold atom physics. 
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Appendix A: Universality within the T-matrix approximation 

To see the universality of the present combined T-matrix theory with LDA, it is convenient 
to write the dimensionless self-energy T,p(iu) n , a s , p,(r), A(r)) = T~ x Yl p {iuj n: r) in the present 
approximation in the form 

POO 

Ep(iu) n ,a s ,fi(r),A(r)) = -T^T ^ / fdq sin 9d0f~ s ' '(iu n ',r) 



xt-.- ~ — - r_y, (Al) 

iuj n+n > - £ q+p (r)T 3 + A(r)Ti 



where p = ^e p /T, q = ^/e q /T, cos6> = p ■ q/(pq), /2(r) = fi(r)/T, A(r) = A(r)/T, 



£ p = £ P /T, and a s = y/2mTa s . In Eq. ( 1A1I) . the dimensionless fermion Matsubara frequency 
u n = u n /T = (2n + l)n no longer has any physical quantity. In the same manner, the 
dimensionless boson Matsubara frequency is simply given by v n = v n jT = 2nir. The 
particle-particle scattering matrix T ~ s (iu n , r) = (27t)~ 2 a/ (2m) 3 TT s q ' (iv n , r) in Eq. flAlj) has 
the form, 



q (.Wn,r) T+ + (iv n ,r) 



lL- r dp - ( nF + ( £i? n,At(r),A(r)) flt + (^ n ,/i(r),A(r)) 
2a s Jo P \Uz-(iu n ,fl(r),A{r)) Ut-{iu n ,fl(r),A(r)) 



T q (iv n ,r) T q + (iis n ,r 

Here, the second term on the right hand side of Eq. ( 1A2|) has been introduced to regularize 



(A2) 



the theory. The correlation function Il| s ' (iu n , fi(r), A(r)) = ((27r) 2 / a/ (2m) 3 T) x IP 8 («V n , r) 
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is given by, for example, 

TL± + (iu n ,Ji(r),A(r)) = jJ2 H P 2 dp sin 6d6-~ ' sA(r) 

s=±l ^° 



X- 



E p+q/2 {r)E p _ q/2 {r) 
£'p+ q /2(r) + sE p - q/2 (r) 



2 



X 



(2mr) 2 + (E p+q/2 (r) + s£ p _ q/2 (r)) 
tanh|4±^^ +S tanh^- /2(r) 



(A3) 



where E p (r) = E p {r)/T = J (p 2 - fi(r)) 2 + A(r) 2 . 

The dimensionless LDA superfluid order parameter A(r) obeys the gap equation, 

1 = --3e f>dp f J- tanh Mi - l) . (A4) 

Thus, A(r) is found to be a function of (a s ,fi(r)). 

Using Eq. (1Alj) . we find that the dimensionless Green's function G p (iuj n ,a s , jl(r)) = 
TG p (icu n ,r) only depends on a s and /2(r), as 

G p -(za; n , a s , /2(r)) = - - - _ . (A5) 

We note that cD n = (2n + l)7r does not involve a physical quantity. In the analytic-continued 
form of Eq. ( )A5j) is given by 



Gp(d) + z5, 5 S , /i(r)) = ~ =— — — = — — , (A6) 

u + iS- £p(r)r 3 + A(a s , /2(r))ri - ^(ui + i5, a s , /i(r)) 

where a) = w/T. As a result, any physical quantity calculated from G p (iu n ,r) can be 
written in the form AT a F(a s , fi(r)), where the coefficient A and the exponent a depend on 
the detailed physical quantity we are considering. F(a s ,jl(r)) is a dimensionless function 
calculated from G, which only depends on (a s ,jl(r)). Thus, when the physical quantity is 
normalized by AT a , it exhibits a universal behavior with respect to (a 8 , fi(r)). In particular, 
in the unitarity limit (where the scattering length a s diverges), the universal behavior is 
dominated by the single parameter /i(r)/T. 
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